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Abstract 

We perform a gauge-transformation on the time-dependent Dirac equation 
describing the evolution of an electron in a heavy-ion collision to remove the 
explicit dependence on the long-range part of the interaction. We solve, in 
an ultra-relativistic limit, the gauged-transformed Dirac equation using light- 
front variables and a light-fronts representation, obtaining non-perturbative 
results for the free pair-creation amplitudes in the collider frame. Our result 
reproduces the result of second-order perturbation theory in the small charge 
limit while non-perturbative effects arise for realistic charges of the ions. 
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I. INTRODUCTION 



Electron-positron pair production from the vacuum in strong electromagnetic fields is a 
fundamental prediction of QED |T|-^|. In heavy-ion collisions at energies near the Coulomb 
barrier, quasi-bound molecular states are formed with binding energies which dive into 
the negative-energy continuum, resulting in a resonance which decays into an electron- 
positron pair In contrast, for ultra-relativistic heavy-ion collisions at peripheral 
impact parameters, the ions execute straight-line trajectories and reside near each other 
for only a very short time. The high charge of the individual ions and the strong Lorentz 
contraction combine to produce fields sufficient for electromagnetic pair production through 
a qualitatively different process. Large cross sections for electro-magnetic pair production 
in these collisions were theoretically predicted || and experimentally observed [|T0| [X7f] . 

In lowest-order perturbation theory, the amplitudes for pair production in heavy-ion colli- 
sions are calculated from two-photon exchange diagrams ||18|| . The quantum field theoretical 
treatment of this process was reduced to a classical- field approach |19[ . Experimental obser- 
vations of free pair production in the energy range around ten GeV per nucleon (in the col- 
lider reference frame) are in agreement with second-order perturbation theory [|IO| , |l1^ ,|14|-|16 



For lower energies of a few GeV per nucleon (collider frame), experimental results for free and 
bound-free pair production show deviations from the predictions obtained from two-photon 



exchange diagrams |11|,|13| . This is likely due to two-center Coulomb effects |j20| , [21|l . In the 
near future, larger ultra-relativistic energies above one hundred GeV per nucleon (collider 
frame) will be available. New non-perturbative effects may become important at colliding- 
beam accelerators such as the Relativistic Heavy-ion Collider (RHIC) at Brookhaven and 
possibly the Large Hadron Collider (LHC) at CERN. These non-perturbative effects are the 
subject of our present work. 

Previous theoretical works on high-energy non-perturbative effects considered unitar- 
ity violation in lowest-order perturbation theory, multiple-pair production, and corrections 
to production cross sections p2Vp2|. An open question, crucial for the beam stability, is 
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whether the non-perturbative effects will enhance or reduce the cross section for bound-free 
production, i.e. production with capture p0| , |32[] . 



In a non-perturbative treatment, starting from the QED Lagrange density operator, the 
Euler-Lagrange equations of motion for the quantum fields are equivalent, under physical 
assumptions, to the one-particle Dirac equation interacting with classical, electromagnetic 
fields J33|. Calculations of probabilities and correlations can then be reduced to solving the 
two-center time-dependent Dirac equation, which describes the dynamics of an electron in 
the classical field of two relativistically moving charges. 

In the ultra- relativistic limit, the ions are practically moving at the speed of light. The 



classical electromagnetic field of a massless charged particle was studied in Ref. p4f . It 
can be described by pure gauge potentials, with different gauges in different regions of 
space-time. The eikonal approximation was then reproduced from an exact solution of a 
quantum-mechanical equation in this field. A similar approach was recently used in target- 
frame calculations of bound-free pair production, where a gauge transformation was used to 
remove the long-range Coulomb effects [ |2"9| , |3"2"| . 

The use of gauge transformations is fundemental to these calculations. The term gauge 
transformation, as used here, is not to be confused with a gauge-symmetry transformation 
in which both the wave function and the fields are transformed so as to keep the equation of 



motion invariant. Here, as well as in Refs. 24, 35,36], for example, a space-time dependent 



phase is used to transform either the wave function or the fields in order to obtain a different 
equation of motion. The connection between the solution of the original problem and the 
gauge-transformed problem depends on the asymptotic (infinite time) behavior of the gauge 
function employed; i.e. on the induced changes to the initial and final states p4| , |35| , |36| . 
Gauge transformations should also be applied with care when used in calculations employing 



truncated basis sets p5|-p7| 



In this work, a gauge transformation is used to solve the two-center Dirac equation 
describing an electron during a relativistic heavy-ion collision. A closed form expression for 
the pair-production amplitudes in the ultra-relativistic limit is found. First, in section II, 



the ultra-relativistic limit for the two-center Dirac equation is obtained and discussed. In 
section III, exact and asymptotic relations between the physical amplitudes and the gauged 
transformed ones are established. In section IV, the light-fronts representation is introduced 
and the foundations are laid for the construction of the exact solution in section V, where the 
perturbative limit is considered as well. The physical contents of our results and an outlook 
for future applications are finally considered in section VI. Details of some derivations are 
given in appendices. 

II. AN ULTRA-RELATIVISTIC LIMIT TO THE TWO-CENTER 
TIME-DEPENDENT DIRAC EQUATION 

Consider a collision between two ions with charges Z A and Z B and velocities (3z and 
—/3z, respectively, moving parallel to each other at an impact parameter of 2b (see Fig. 
|l]). An external- field approach to the influence of these ions on the vacuum is appropri- 
ate for peripheral impact parameters (i.e. no nuclear interactions), heavy ions, and ultra- 
relativistic energies, when to a very good approximation, the ions continue intact on their 
parallel, straight-line trajectories. The two-center Dirac equation for an electron in the 
time-dependent external field of these ions is given by: 

i^imt)) = [H + H A {t) + H B (tj\ |$(f,t)), (1) 

where |$(r, t)) is the Dirac spinor wave function of the electron, H is the free Dirac Hamil- 
tonian and HA{t) and H B {t) are each the interaction with one ion, 

H = -ia ■ V + 7°, (2) 

H A (t) ee (I 4 - (3a z ) ~ ZAa (3) 

V(r± - b?/l 2 + - pty 

H B (t) ee (I 4 + (3a z ) . ~ ZBa (4) 

^{f ± + b)yj 2 + {z + (3ty 

We are working in the collider frame, using natural units (c = 1, m e = 1, and % = 1), 
and applying the conventional notation; (3 = v/c, 7 = 1/Vl _ P 2 - ot is the fine-structure 



constant, at and 7 M are Dirac matrices in the Dirac representation, as in Ref. [p]]; a are the 
Pauli matrices; and I 2 , 2 , I4, and 4 are the 2-dimensional and 4-dimensional unit and zero 
matrices. 

We would like to consider the ultra-relativistic limit in which 



0->l, 7>&,^±- 



(5) 



Eq. (0) does not simplify in this limit in a straight-forward way because, for any given 
time, the long-range behavior of the interaction terms H A and H B is independent of 7 as 
z — > ±00 (see Fig. |^). A simple ultra-relativistic limit can be obtained by first applying a 
gauge transformation in order to remove this long-range tail of the interaction. The gauge- 
transformed wave function \^{r, t)) is defined by 



\*(?,t)) = U(z,t)\$(r,t)), 



(6) 



U(z, t) = exp \%Z A a In -7^ - z) + ^b 2 + 7 2 (t - zf 



+ iZ B a In 



+j(t + z) + ^b 2 + 1 2 (t + z) 2 



The Dirac equation for |\l/(r, t)) is obtained from Eq. ([[]), 

i^(f,t)) = [H + W A (t) + W B {t)] mr,t)) 



where the new interaction terms are 

W A (t) = (l 4 -(3a z 

-(I 4 - a z 

W B (t) = (l 4 + (3a z 

-(I 4 + a z 



-Z A a 



(r± - b) 2 h 2 + (z- (3t) 2 
-Z A a 



fe2/ 7 2 + ( Z -t) 2 

-Z B a 

(r-l + 6) 2 / 7 2 + (z + (3t) 2 
-Z B a 



(7) 



(8) 



(9) 



(10) 



l b 2 / 1 2 + (z + t) 2 

Figure |3| demonstrates the short-range character of this gauge-transformed interaction. Sim- 
ilar gauges have been used in Refs. |38|]24|] which reduce in the limit /? — > 1 to Eq. (0). Unlike 



Eq. (JT]), the gauge-transformed equation, Eq. (||), has a simple ultra-relativistic limit fl3~2"| , |34" . 
In the limits of Eq. (|^) (large 7, small r±, and small impact parameter &,), Wa ( Wb) has a 
sharp, delta-function dependence on t — z (t + z) p2| , pl| , p9| , (see appendix [A]), i.e. 



W A -> (I 4 - a 2 )Z4C*5(t - 2) In 



-> (I4 + a z )Z B a6(t + *) In 
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(rl + &) 2 
6 2 



(11) 
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Consider the physical nature of this limit. A 5 function over time alone would indicate 
a sudden interaction of the ions with the vacuum. In the gauge-transformed equation, 
Eq. (§), with the interactions of Eqs. ([□]) and (0), as they move, the ions are continuously 
interacting with the vacuum. Naturally, this interaction is singular on the trajectories of the 
ions, as it was before the ultra-relativistic limit has been taken; but an additional singularity 
is induced in the ultra-relativistic limit by the extreme Lorentz contraction of the field. In 
this limit, the interaction is infinite on the two planes perpendicular to the ions trajectories, 
and vanishes elsewhere. In the following, we calculate pair production amplitudes using 
Eq. (H) with the interactions in Eqs. (|i~T|) and (|12"D . The region of large r± is not properly 
accounted for in this treatment, but contributions from this region to pair production are 
assumed to be small. 

The interactions in Eqs. (|TT| ) and fll2|) have zero range in the longitudinal direction and a 
logarithmic behavior in the transverse direction, similar to the potential of a line of charge. 
In the limit (3 — > 1, the two ions are moving at the speed of light and thus the interaction 
planes described above coincide with the light fronts, given by z = ±t (see Fig. Finally, 
we note that (I4 ± a z )/2 are orthonormal projection operators. The 4-Dirac spinor wave 
function of the electron can be decomposed into two orthogonal components, 



l* + (^*)> = ^(U + & z )\W,t)) 

|*_(f,t)) = i(l4-& z )|*(f,f)). 



(13) 
(14) 



Each ion interacts directly only with one of these components; Za with |\I/_(r, t)) and Zb 
with \^ + (r,t)}. 



III. ASYMPTOTIC SOLUTIONS AND TRANSITION AMPLITUDES 



In scattering theory, characterized by free initial and final states, a complete solution is 
generally given by the set of asymptotic transition amplitudes between plane waves, (the 
S-matrix). In this section, we define the transition amplitudes, S% , for the electronic 
spinor wave function |$), i.e. for Eq. (|I|), and the transition amplitudes, for the gauge 
transformed wave function i.e. for Eq. (|8|), using, as usual, initial and final plane wave 
states. Use of an initial condition of a single plane wave for Eq. ([!]) is somewhat subtle since 
the theory is never free due to long-range Coulomb effects O. This is an important point 
when integrating the amplitudes to obtain predictions for physical observables. One should 
properly account for actual initial and final states for a given experiment. 

A complete set of solutions of the free Dirac equation is given by the Dirac plane waves; 
{|x p (r, t))}. Each plane wave is characterized by three continuous and two discrete quantum 
numbers; namely, the three components of the momentum, p, the sign of the energy, (we 
use A p = for positive energy electrons, A p = 1 for negative energy electrons,) and the spin, 
i\s p ) = |+) for spin up and \s p ) = |— ) for spin down,) p = {p,X p ,s p }. The plane waves 
satisfy, 



H \x P (r)) = E p \xp{r)), 
\x P (r,t)) = exp(-iEpt)\x P (r)} , 



E„ 



{-iy p e p , e p = Jp 



(15) 
(16) 
(17) 



where Eq. ( |17D is the condition for being on the energy shell. An explicit form is given, for 
example, by [[|, and in our notation by 



\Xp\r)) = i \ exp^r • p) 



2e p (l + Bp) 



u 


-0 




v 1 * 


2 J 


\ 



l+Sp) \s p ) 
a ■ p \s p ) 



= exp(zr • p)\u p ) . 



(18) 
(19) 



We define the solution \<f>^\r,t)) of Eq. ([!]) by the initial condition, 

lim \^\f,U)) = \Xj(r,ti)). (20) 



ti— > — oo 



The asymptotic transition amplitude S^' is then given by 

S?>= lim ( Xk (f,t f )\^\f,t f )), (21) 

t/-+oo 

where the bra-ket stands, as usual in non-relativistic and non-covariant notation, for inte- 
gration over all space fat a given time. Likewise, \ip^(r,t)) is defined as the solution of 
Eq. (^|) with the initial condition, 

lim \^\f,U)) = \ Xj (r,U)), (22) 
and the asymptotic transition amplitude is given by 

A<£>= ]im{ Xk (f,t f )\^ j) (r,t f )). (23) 



The initial condition for Eq. (H), Eq. (|20[) , and the initial condition for Eq. (||), Eq. (p2|). 
correspond to different initial physical states, as they are not related by the gauge trans- 
formation in Eq. ([F|). A similiar comment is true for the final states used in defining the 
amplitudes in Eqs. fl2T|) and fl23|) . In general, S k and A^' are completely different ampli- 
tudes. They are related to each other by the gauge transformation of Eq. ([?D in the following 
way, 



CO = 

xA® (XqftUWMlXjftU)), (24) 



Si 3) = E E (Xk(r, t f )\U\z, t f )\ Xp (r, t f )) 
p q 

i( ( 

p 



where J2 P stands for integration and summation over all the quantum numbers, p = 
{p, X p , s p }. This relation is based on the completeness of the plane- waves basis set and 



should be questioned if applied with a truncated basis calculation |35|-|37|1. 

A relation like Eq. fl2~3| ) holds between any two amplitudes which are related by a gauge 
transformation. If the gauge transformation U had been otherwise defined so that it would 
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become unity for asymptotic times, (i.e. U — > 1 for t% — > — oo and i/ — > oo), the orthonor- 
mality of the plane- waves would eliminate the double sum in Eq. (123), and the asymptotic 



transition amplitudes and would be identical. Gauge transformations which share 
this property have been described as exhibiting asymptotic gauge invariance [p^j . Like- 
wise, if the gauge transformation U would be independent of r for asymptotic times, (i.e. 
U — > exp(iCj) for ti — > — oo, and C/ — > exp(iCf) for — > +00, where C!j and C/ are real 
numbers), then the double sum in Eq. (|24l) would again be eliminated and the matrix ele- 
ments of U in Eq. fl24|) would contribute only a single constant phase at asymptotic times. 
In this case, the asymptotic transition probabilities derived from Sj? and would be 



identical |2j 

The specific gauge transformation used here, Eq. (|7|), is not independent of space at 
asymptotic times, as is shown in appendix |Bj, and, as a result, does not exhibit asymptotic 
gauge invariance for the amplitudes or the probabilities derived from them. However, this 
gauge does have an additional property which relates Sj? and in a way more useful 



than Eq. fl23|). It is shown in appendix [TJ that in the special case of symmetric collisions, 
with Za = Zb, can be expressed as a series expansion in powers of 1/tf and l/\ti\ whose 



zero-order term is A k 



IV. THE SHARP DIRAC EQUATION IN THE LIGHT-FRONTS 

REPRESENTATION 

In this section, the sharp Dirac equation, Eq. (H), with the limiting form of the interaction 
in Eqs. (|TTD, and (0), will be further simplified by changing into light-front variables and 
by introducing a new representation for the Dirac spinors, the light-fronts representation. 
This is an appropriate choice of variables and representation, since, in the ultra-relativistic 
limit of Eq. (|5|), the interactions are confined to the light fronts. 
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A. Definitions and notations 



In terms of light-front variables, space-time and energ; 
4- vectors (r_j_,r + ,r_) and (p±,p + ,p_), where 

r± = (t±z)/2 
p± = E p ± p z 



P+P- 



i+Pi 



;y-momentum are described by the 



(25) 
(26) 
(27) 



The sign and absolute value of {p + +p_)/2 are \ and e p , respectively. Equation ( |2"7| ) like 
Eq. (|I7D defines the energy-shell. These variables were often used previously for quantization 



on one of the two light fronts, r + = or t_ = ||40|1 . For the problem considered here, it is 
useful to keep the symmetry between r+ and r_. 

The projection operators (I4 ± a z )/2 acquire a simple form and the interaction is diago- 
nalized by introducing the light-fronts representation for the Dirac matrices, 



7ught- fronts — ^TDirac"^- ) 

A = 

V2 



Aa^A t 



(l 



2 2 
2 -l 2 



A 



A 



-(I4 + " 2 



At 



At 



( h0 2 \ 



y 2 2 y 
' o 2 2 » 



v 2 I 2 j 



Ad 1 A f 



UJ O5 



(28) 
(29) 

(30) 

(31) 

(32) 

(33) 
(34) 
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With this notation, the gauge-transformed two-center Dirac equation in the sharp ultra- 
relativistic limit in the light-fronts representation is 

t ■> ■ ^\ { xfn-.\nf*A h~\ f \n.\\ 



id T+ \G + ) 
id T _\G^) 



5(r + )B(r ± ) h 



G + ) 



(35) 



where \G+) and are the upper and lower bi-spinor components of the Dirac wave 

function in the light-fronts representation 

\ 

= A|¥), (36) 



G+) 
G-) 



and 



h = l 2 -iuJ- p ± , 
A(tj_i b) = Zaoc In 



B(f±, b) = Zboi In 



(gi. - b) 2 
b 2 

{r ± + b) 2 
b 2 



(37) 
(38) 

(39) 



The upper and lower bi-spinors are coupled by the free Hamiltonian. Each interacts directly 
with the external field of one ion and feels the field of the other ion through its coupling to 
the other bi-spinor. 

Equation fl3"5|) has no discontinuities in the transverse direction. It is therefore useful 
to Fourier transform its solution with respect to r±. Two mixed bi-spinors wave-functions, 
\g±(q±; r + , r_)), are then defined by 



\G±(r±,T+,T-)) 



,jr ± -q A 



\g±{q±;T + ,T_)). 



(40) 



|^ + ) and 1^-), like \G+) and are coupled by the free Hamiltonian. 



B. Free Dirac equation off the light fronts 



Off the light fronts, i.e. for r + ^ and r_ ^ 0, the wave function satisfies the free 
Dirac equation and Eq. (j3~5"D reduces to two coupled equations for the mixed bi-spinors 
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\g±(qr,T + ,T..)). 



d 

% -Q^r \9+) = (h - iu ■ <t±)\g-), 
d 

OT- 



(41) 
(42) 



As usual, the second-order equations decouple 

d 2 



dr+dr- 



g ± ) = -(l + q 2 ± )\g ± ), 



(43) 



where use was made of 



(I 2 - iuj ■ q±)(h + id • q±) = (1 + q]_)h- 



(44) 



A solution to Eqs. (0, H) 

is given, for example, by the plane waves of Eq. fll8|) which 
in the light-fronts representation are given by 



( \F%) ^ 



\FQ = J dq ± e^\f ± (q ± ;T + ,r_)), 

\fUqx; t+,t_)> = S(q ± -p ± ) e -^++^-)|Il). 

The bi-spinors, |r±), (the upper and lower parts of A\u p )), 

(2tt)- 3 / 2 



(45) 

(46) 
(47) 



x 



2je p {l+e p ) 



h (l + (-l) Xp p±) ±id-p ± ] (±a z ) x *\s p ), 



(48) 



satisfy the simple relation 



= I 2 -iu-p ± 

P+ + ' 



(49) 



These plane waves solve Eq. (pop off the light fronts in the limits t — > ±oo. They do not 
solve it for finite i, when pj_ is no longer a good quantum number, as the singular interaction 
with the ions makes the wave function discontinuous at the light fronts. 
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C. The discontinuity across the light fronts 



It is standard procedure in wave-mechanics to form piece-wise solutions by satisfying 
continuity relations at the boundaries between free regions. It was shown in Refs. |]32|j34l| 
that a 5-function singular interaction at a light front results in a discontinuity in the electron 
wave function which is given by a space-dependent phase shift. The proof is reviewed in 
appendix 0, where it is shown that for our case of Eq. fl35l), the discontinuity is 

\G + (r + = 0+)> = e- iB ^\G + (r + = (T)>, (50) 
|GL(r_ = 0+)) = e- M ^|GL(r_ = 0-)). (51) 

These phase shifts are derived from Eq. fl35|) in general for any A(f±,b) and B(f±,b), i.e. 
any functional dependence on the perpendicular coordinate. Here, A(f±, b) and B(f±, b) are 
given by Eqs. (|38| , p9|) . 

Due to this space dependent phase-shift, the transverse momentum is not conserved and 
the Fourier components of Eq. (|40"|) are mixed when the singularities at the light fronts are 
crossed, 

\g+{q±] r+ = 0+)) = J dp±Qz B (p±_ - q±, -b) 

x\g+(pr,T + = Q-)), (52) 

\g-(q±; r_ = o + )) = J dp±Q ZA (p± - q±, b) 

xb_(p ± ;r_ = 0-)), (53) 



where 



(r± ~ b) 
b 2 



2 laZ 



(54) 



Note that here k and b are two-dimensional vectors in the (x, y) plane. The continuity 
is recovered in the limit Z — » 0, as Qz(^,b) — > 8(k). The distribution b) in general 

diverges. This divergence is an artifact of applying the sharp limit for the gauge-transformed 
interaction, Eqs. flTTD and (0), for large r±, i.e. outside its range of validity. The properties 
of the distribution Qz{k, b) for finite charge are considered in appendix 0. 
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V. A PIECE- WISE SOLUTION TO THE SHARP DIRAC EQUATION 

In this section, the formalism that was introduced in section IV is used to obtain the 
transition amplitudes between asymptotic plane waves, A^, defined in section III. 

The singular interaction on the planes perpendicular to the trajectories of the ions, cut 
space-time along the light fronts into four regions, as is shown in Fig. [5[ A piece-wise 
solution is defined off the light fronts by \g±(q±; t+, r_)) = \g±{q±; t+, t_)), where (i)= I for 
r + < and r_ < 0, (i)= II for r + > and r_ < 0, (i)= III for r + < and r_ > 0, and 
(i)= IV for r + > and r_ > 0. In each region, the wave function is continuous and solves 
the local free Dirac equation. At any time, except for t — > ±oo, the wave function extends 
in space through three (or two, at t = 0) of these regions. The solution presented here is 
not complete in the sense that it does not include the solution on the light fronts; r + = 
and r_ = are excluded. The physics on the light fronts may contribute to bound-free pair 
production. Thus, our present work is limited to free pair production. 



A. Initial condition and intermediate states 



Consider the initial condition, Eq. fl2"2"|), of a single plane wave with the quantum numbers 
j = {j, Xj, Sj}, or, using light-front variables, j = {j±,j+,j~, Sj}, with the constraint j+j- = 
1 + The continuity off the light fronts gives the solution in region I, 



\gi(q x )) = 5(j ± - gi)e-^+ + ^-)|ri>, 



(55) 



where the bi-spinors \T±) are defined as in Eq. (|48|) . 

The solution in regions II and III is obtained by first applying Eq. ([52]) for the disconti- 
nuity across r + = and Eq. (|53D for the discontinuity across r_ = and then solving the 



coupled equations ( |41| , [42]) inside each of the intermediate space-time regions. We obtain in 
region II 



\g+(q±)} = exp 



' i + gT 
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xQz B (j±-q±,-b) |L J + ), 
- iu 
3+ 



i«E«o> = ( | ff ; (9 - ±)) , (56) 



and in region III, 



9 IU (P±)) = exp 



3- 

xQz A 0±-px,b) |ri), 



|0?(£l)> = ^-^'^ j |^(gl)). (57) 

It is now apparent why the Fourier transform with respect to r± and the definition of 
\g±(q±; r + , r_)) in Eq. PD| ) were needed. The simple discontinuity condition fl5"0|) at r + = 
applies only to The other bi-spinor \G-) is influenced indirectly by the field at r + = 

through its coupling to \G + ). Likewise, at r_ = the simple discontinuity condition (51) 



for induces a non-trivial change in \G + ). The coupling between \G + ) and in free 
space on either sides of the singular interaction is best described by Eqs. (41,42) for their 



Fourier components with respect to f±. Thus, while the discontinuity conditions (|5^J53| ) for 



\g±) seem more complicated than the discontinuity conditions (|50| , |51[) for \G±), using \g±) 
allows for a simple derivation of the complete spinor wave function in regions II and III. 

It is a well known fact that two ions are needed in order to create an electron-positron 
pair. This can also be seen here. In the presence of ion B alone, for example, Eqs. (|56D give 
the solution for r + > 0, including the asymptotic solution at tf — » oo. Projection on a plane 
wave then gives a conservation law for the positive light-front momentum. Likewise, in the 
presence of ion A alone the negative light front momentum is conserved. 

A^(Z A = 0)^5(k + -j + ), (58) 
A<l\z B = 0)<x6(k--j-). (59) 

A direct result from k + = j + or k- = j- is that the sign of the energy of the electron is the 
same before and after the collision. Thus, our formalism satisfies the known result that the 
passage of a uniformly moving charge does not induce a transition changing the sign of the 
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energy. The presence of both ions defines a new region of space-time, region IV (r± > 0) 
which is the space between the ions (— t < z < t) after the collision (t > 0), i.e. when the 
ions are already moving apart. It is shown below that the sign of the energy can change and 
pairs may be created in transitions from the initial state in region I (i, — > — oo) to the final 
state in region IV iff — > oo). 

The solution of the free Dirac equation in region IV is complicated by the non-trivial 
boundary conditions on the light fronts. Applying Eq. ( |52f) again for the discontinuity across 
r + and Eq. (|53|) for the discontinuity across r_, we cross from regions II and III into region 
IV to obtain on the hyper- surfaces adjacent to the light fronts, 



<I<1 cxp — /r_ 

x Qz A (q± - k±,b) 

x Qz B (j±-q±,-b) 
( I 2 + iuj ■ q± 



3+ 



3+ 



(60) 



j dp± exp 


Wf)] 







x Qz B (p± - kj_ - b) 



x Qz A {3L-p±,b) 
( I 2 - iuo ■ p±\ 



x 



J- 



ri). 



(61) 



Instead of solving now for \g±) at any r± > 0, the transition amplitudes are obtained in the 
next subsection by defining the transition current and by applying Gauss' theorem for this 
current. 

B. The amplitudes 



The transition amplitudes were defined in Eq. 



4 J) = lim / drxl(r,t f ) ^ } {f,t f 



(J), 



(62) 
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The integrand is a component of a 4-vector current density, which is a conserved quantity 
(see appendix 0). This transition current P5[ , (J*^'^, Jq^), is defined by 

■t 3) = xl ^ 

J^>= (63) 

An equivalent form for the transition current in terms of light-fronts representation wave- 
functions includes 

jgA = jM) ± jfrj) = 2F^G i £\ (64) 

It is now possible to use Gauss' theorem on the hyper-surface of the inner border of 
region IV to show that (see appendix [F|) 

Af = lim / dr ± f°° dz 4 kJ \r,t f 



tf^CO 



— OO 



2 J df ± J™ dr_ 4 M (f ±) r + = + , r_ 



-2 J df x J q+ dr + J ( y\f x , t+, t- = 0+). (65) 

The transition currents J±'^ are calculated from the results of the last subsection by using 
Eqs. (|gg|443,||). 



J± ,3 \f±,r + ,T_) = 2 J dp±J dl±exp[if± ■ (l±-p±)} 

x (fi(p±; r +l r_M(l ±] t+, r_)>. (66) 

Integrating over f± and using the explicit expression (^) for the plane waves, 
Af = rdr^- k +(V k + \g^(k x ;r + = 0+ T-)> 

roc _, 

-2 / dT + e iT + k -(r k _\g™(k r ,T+,T„ = 0+)). (67) 



The amplitudes are finally obtained by substituting Eqs. (60,61) and integrating over t±. 
The integration over r± would have given a 5-function conservation law for the light-front 
momenta, had it been on the complete line — oo < r± < oo. Instead, the integrals on the 
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half lines < t± < oo are regulated in the usual way with an infinitesimal small constant, 



V 0- 



— 'i 

drexp(irK) = . (68) 

0+ K + 17] V ' 



The transition amplitudes corresponding to the exact solution of the sharp Dirac equation 
off the light fronts are 

a? = x 

in 

I J dp ± Qz B (p±_ - k±, -b) Qz A (l± - P±, b) 

(T k + \I 2 - i£ ■ p ± \Tl) 
j-k+ - (1 + p\) + ir](-l) x J 

- / dq ± Q Za {q± - k±,b) Q Zg (j± - q±, -b) 

x mh+j^n) \ 

i+k. - (i + ?i) + i,,(-i) x > / • 1 ' 

where the infinitesimal small, positive constant, rj, can be omitted for pair-production am- 
plitudes corresponding to Ej < and Ef. > 0, i.e. j±k T < 0. This is an interesting result. In 
the ultra-relativistic limit of Eq. (|5[), the asymptotic time evolution of the gauge-transformed 
electron wave function, \ty(r,t)), is exactly given by these amplitudes. An exponential, non- 
perturbative dependence on the coupling constant aZ appears here as non-trivial phases 
in the integral representation of the distributions Qz A (i^±,b) and Qz b (k±, —b) which were 
defined in Eq. (0). The two terms in Eq. fl69| ) correspond to two different time-orderings of 
the interaction with the ions. In the next subsection, they are shown to reduce in the small- 
coupling perturbative limit to the well known two photon exchange diagrams as depicted in 
Fig. [7|. For vanishing charg could be expected, 

Af(Z A = 0, Z B = 0) = 8{k - j)6 Xk>Xj 6 Sk>Sj . (70) 



C. The perturbative limit 



The small-charge perturbative-limit of the pair-production amplitude was calculated in 
Ref. [18). To leading order in aZ (second order), the amplitude is given by a sum over two 
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diagrams, where each diagram describes a two-photon exchange process. The second-order 
perturbation-theory result, S% , for the transition amplitude between an initial negative- 
energy state j = {j,Xj = 1, Sj} and a final positive-energy state k = {k,\k = 0,Sk}, is 
given by Eqs. (24-32) of Ref. |IBfl . In the ultra-relativistic limit, (3 — > 1 and large 7, the 



perturbative result reduces to 



x 



$1 (j) = J d P_L exp[-ib ■ (2pj_ - j± - k±)] 

18 (aZ A )(aZ B ) 
x ^ — — 

(rl - k±) 2 (p_L - j±) 2 

(T k + \i 2 -t£-p ± \r j _) 

3-k + -(l+p 2 ± ) 
J dq_i_exp[ib- (2q± - j± - k±)] 

i8 (aZ A )(aZ B ) 
(q± ~ k±) 2 (q±-]±) 2 
(V k _\I 2 + id-q x \V j + ) 



x 



X 



j + k_-(l + ql) ■ (71) 

A transformation to the light-fronts representation was used here to obtain bi-spinor bra-kets 
from the 4-spinor bra-kets of Ref. QT8 1 . For example, using At A = I 4 , 

(u k \ (h - ocz) (& ■ Pi. + 7°) {h + & z ) \uj) 

= (T k + \I 2 -tO-p ± \rl). (72) 



It is interesting to compare the perturbative result of Eq. ( |7T| ) to our non-perturbative 
result of Eq. fl69"|). In the small-charge limit of aZ — > 0, after proper regularization, the 
leading-order perturbative limit for Qz from appendix |D| can be used, 

r> ^ A r/^ ^Z 1 r . r 

Qz(K,b) —*■ o{k) -expzO-zc. (73) 

7T K 

Direct substitution shows that in this limit the non-perturbative result of Eq. (|69| ) exactly 
reproduces the perturbative result of Eq. (|7ip . 
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VI. CONCLUSIONS AND OUTLOOK 



We have used a gauge transformation to obtain a useful ultra-relativistic limit for the two- 
center Dirac equation, which allows for an exact solution off the planes perpendicular to the 
ions' trajectories, i.e. off the light fronts. In general, the amplitudes of the gauge-transformed 
Dirac equation are related to the amplitudes of the original equation in a non-trivial way 
due to long-range Coulomb effects. For symmetric collisions, and calculations in the collider 
frame, some of these long-range effects cancel. The two different amplitudes are then related 
by a series expansion, and, to leading order, they are equal. 

The amplitudes were calculated here in the ultra-relativistic limit, assuming 7 to be 
large. No assumption was made on the value of the charge times the fine-structure constant 
Za. When taking the limit of small Za, we are able to show a complete agreement with the 
ultra-relativistic limit of the expression obtained from standard second-order perturbation 
theory [|L8[| . In second-order perturbation theory, pair production is described as a two 
photon exchange process in which each ion exchanges one photon with a negative-energy 
electron. The negative energy electron is kicked off its energy shell by the first interaction 
and then kicked back to the energy shell by the second ion, but with a positive energy. 
The two diagrams that contribute to the amplitude differ in the time order of these photon 
exchanges, or 'kicks'. Our result, which provides a very similar physical picture of pair 
production as a 'two-kicks' process, is obtained in the ultra-relativistic limit within a rather 
different, and completely non-perturbative approach. 

In our work, the electromagnetic fields of the ions are confined to the light fronts by 
the extreme Lorentz contraction and by a choice of a particular gauge designed to remove 
the long-range Coulomb effects. In this gauge, as the velocity of the ions approaches the 
velocity of light, each ion carries with it, perpendicular to its trajectory, a wall of singular 
electromagnetic interaction. An initial plane wave in the space between the approaching 
ions acquires a space-dependent phase shift as it is swept by this singular-interaction wall. 
A single plane wave between the ions gives a distribution of local plane waves in the space 
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behind each ion. Had there been only one ion, no transition would be allowed between 
the negative-energy continuum and the positive-energy continuum, i.e. no pairs could be 
produced. Pairs are produced because, as the ions move past each other, the two phase- 
shift planes collide. After the peripheral collision, as the ions move apart, the solution in the 
space between them is determined by the non-trivial boundary conditions at the light fronts. 
The main result of this work, the exact integral representation for the free pair creation 
amplitudes of Eq. (|69D , is finally obtained in this framework by calculating the transition 
currents flowing from the light fronts into the space between the separating ions. The two 
terms correspond to the two time orderings of the interaction of the two phase-shift walls 
with the electronic wave function. In the perturbative limit of a small coupling constant, 
the effect of the singular field perpendicular to each ion reduces to a single photon exchange. 
For finite charges of the ions, the perturbative linear dependence of the amplitudes on each 
charge is replaced by non-perturbative, nontrivial phases in our integral representation. 

Numerical evaluation of the non-perturbative effects, differential cross sections, and ap- 
plications to multiple pair production will be considered in future work. 
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APPENDIX A: THE SHARP, ULTRA-RELATIVISTIC LIMIT 

In this appendix, we will outline the derivation of the 5-function limit of the electro- 
magnetic interactions, Wa and Wb, given in ( |TTJ ) and (|E^) beginning with their definitions 
given in ([]) and (PUD, respectively. (The same limit has been previously obtained; see Refs. 

3|@lg3). 
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In the limit of extreme ultra-relativistic collisions, one may neglect terms in the interac- 
tion proportional to 7~ 2 . It is then possible to first set (3 — > 1 and then use 



b 



2 v-1/3 / 2 \-V2 



■^fHr)u(t). (A1) 



6 2 



It is easy to verify that in the limit 7 S> a, 6, 



and that in the same limit 



dr V (r;a,b)^\nl-\, (A2) 



a 2 - b 2 

T](t ^0;a,6)cx >0. (A3) 

7 2 



APPENDIX B: THE GAUGE TRANSFORMED AMPLITUDES 

In section III, the transition amplitudes and were defined and an exact relation, 
Eq. (P^), was established between them. However, this relation is not always useful as it 
involves infinite integrals and sums over the complete plane- wave basis set. In this appendix, 
a series expansion in inverse powers of the asymptotic time will be shown to relate Sjf and 
in a simpler way for symmetric collisions and calculations in the collider frame. 

Two evolution operators can be defined for |<3>) and |\&), from Eq. (^) and respectively, 

|$(f,t / )) = V(t / ,t l )|$(f,t l )) (Bl) 
\*(?,t f ))=U(tf,ti)\W,U))- (B2) 

They are related by the gauge transformation of Eq. (|7|) 

U(t f ,ti) = U(z,t f )V(t f ,t t )U\z,t t ). (B3) 

The amplitudes are given by 
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Sjp = lim (Xhfttf)\V(tf,U)\Xi(?,ti)), (B4) 

lim (Xfc(r,t/)|W(t/,ti)|^(r,ti)>. (B5) 
*/ 

A direct substitution gives 



t f — >0O 



t* — too 



S { k J) = hmJ X k(r,t f )\U\z,t f )U(t f ,U)U(z,U)\ Xj (f,U)). 

(B6) 



The asymptotic expressions, 

U\z,t f 



^ [27(1^1-^/6]^°' 
. t^oo [2j(t f - z)/&]^ a 



[27(t/ + z)/&] iZ * a ' 
reduce, for symmetric collisions, (Z = = Zg), to power series in z/t 

U(z,U) ~ 1 + i2Za-^- + ... 

\h I 

C/ f (^t/) « 1 -i2Za— + .... 
Substituting these power series in Eq. (|B6|) and integrating term by term one gets, 

Ou ~ flu 



(B7) 
(B8) 



(B9) 
(BIO) 



+i2Za lim (Xk(P,tf)\U(t f ,ti)r-r\xj(r,ti)) 



-z2Za lim ( Xfc (f,t / )|fw(t / ,t i )|x,(f,t i )) + - 

t/ — «XJ 



Using completeness, Eq. ( |B11| ) can equivalently be written as 



(Bll) 



+i2Za J2 



4° lim 

I j — > OO £ - 



- Ap 5 lim (x* (f , t/) I f |xj (r, t») > 



t f — >oo 



+ 



(B12) 



It is clear that and are in general different. For non-symmetric collisions, the 
relation between them involves, for example, a highly oscillatory, z-dependent phase which 
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explicitly depends on 7. However, this phase cancels for the interesting case of symmetric 
collisions. The first-order corrections to ~ A^' decrease linearly with time and have 
a functional nature of polarization effects. Higher-order corrections are higher moments of 
z/t. 

APPENDIX C: THE DISCONTINUITY AT THE LIGHT FRONTS 

In this appendix, the discontinuities of the spinor wave function at the light fronts (at 
r + = and at r_ = 0, excluding only r + = r_ =0) are deduced from Eq. fl35l). Previous 
derivations of the discontinuity of a wave function due to an ultra-relativistic charge are 
reviewed. 

At one light front (t+ = 0, r_ ^ 0), Eq. (|35|) for \G+) reads, 

id T+ \G+) = h \G-)+B(r x )6(T + )\G + ) . (CI) 

The 5-function singularity renders \G + ) discontinuous at r + = 0, as can be seen by inte- 
grating both hand sides of Eq. (|Cl| ) with respect to r + from — e to e and taking the limit 

\G + (r + = + )> ^ \G + (r + = 0-)). (C2) 
An auxiliary bi-spinor can be defined by a piece-wise gauge transformation, 

\G + )=exp[iB(r ± )9(r + )]\G + ). (C3) 

Direct substitution gives, 

id T+ \G + ) = e W [iB{r ± )6(T + )]h \G„) (C4) 

The auxiliary bi-spinor is continuous at r + = 0, as can be seen by operating on both sides 
of Eq. ( |C4| ) with lim e ^ fl e dr + , obtaining 

\G+(r + = + )) = |G + (r + = 0-)). (C5) 
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The continuity of \G+) at r + = (r_ 7^ 0), implies the discontinuity of Eq. ( |50D for \G + ). 
Likewise, a continuity of 

\G-)=exp[iA(r ± )9(T-)]\G-) (C6) 



at t_ = 0, (r + 7^ 0,) implies the discontinuity of Eq. ( pT| ) for 

This Heavyside step-function, space-dependent, phase discontinuity was previously ob- 
tained in Ref. f32| . In earlier work [34j], a gauge transformation was used to establish the 



fact that the electromagnetic field of a charge which is moving at the speed of light can 
be equivalently given by gauge potentials with a 5-function singularity at the light front, 
or by gauge potentials with only a step-function discontinuity there. The wave function 
of a particle interacting with this field is discontinuous or continuous, depending on the 
gauge choice. We choose to work with such a gauge that would give a sharp interaction 
and a discontinuous spinor wave function, yet we have used here other gauges to find the 
discontinuities in an explicit form. 



APPENDIX D: THE MOMENTUM-TRANSFER DISTRIBUTION 

When the singularities at the light fronts are crossed, the transverse momentum changes. 
The distribution for this momentum change is given in section IV by Eq. (p4|), 

-iaZ 

(Dl) 

As mentioned in section IV, divergence of Qz{H,b) is an artifact of the ultra- relativistic 
approximation used in Eq. (|TT|) and (|i"2"|). The integral over this distribution converges and 
is normalized to 1 

JdZQ z {K,b) = 1. (D2) 

For a vanishing charge, Qz=o{H,b) = 5(H), but for finite charge Z / 0, this distribution 
diverges both for vanishing and finite momentum transfer k. In this appendix, we show that 
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Qz{Z,b) = 



(2tt)* 



dr\ e" 



(fj_ - b) 5 
b 2 



for H 7^ in the perturbative limit, aZ <C 1, and after proper regularization, the leading 
order correction to the 5-function is given by Eq. (|73|), i.e. 

%(y.Z 1 

Qz(H,b) —y 5(H) ——exp(ib-H). (D3) 



7T K 



Integrating first over the angular variable, 

Qz(H * 0, b) = &2eXP ( ^^ } r dsMsbnjs 1 -^, (D4) 

2,71 Jo 

where b = \b\, k — \H\, and Jo is the Bessel function. The integral over s diverges, but can 
be regulated for finite k in the limit of aZ <C 1. 



Using, for example, Eq. (6.631.1) in Ref. [42] and Eq. (13.5.1) in Ref. |E| one has 



1 f°° 2 

lim— / dse es Jo[sbK]s 
e-+o 2ix Jo 



\-%2aZ 



r r(l-iaZ) / -(bK) 2 \ 
1 (ibK\~ i2aZ 

= !"Ste exph(6,c)2/4e| (irj 

-!?V)-V-^ fe)*', (D5) 
7r T(+iaZ) \ 2 J 

where \F\ is the confluent hypergeometric function. Eq. ([73]) is now obtained by taking 
aZ — > and using 5(6/5) = lim e ^ ^ exp[— (6«;) 2 /4e]. Note that the limit e — ► can only 
be taken after taking the perturbative limit aZ — ► 0. 



APPENDIX E: THE TRANSITION CURRENT 

As we are unaware of an appropriate reference, we prove in this appendix that the 
transition 4-current density defined in Eq. (|63|) is conserved. In fact, any two solutions of 
the free Dirac equation can be used to define a conserved current in a similar way. This 
proof is very similar to the one found in textbooks proving the probability current to be 
conserved [0. 

Both Xk and ip^ solve in region IV the free Dirac equation in the Dirac representation 
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d - 

i —if)V>(r,t)= -id -V + 7°JV (i) (^*) 



■9 V 



-ia-V + 7° Xk(r,t) . 



(El) 
(E2) 



Multiplying Eq. 



from the left by the adjoint of Xk gives 



dt 



U) _i_ -vt-Ai/iO') 



(E3) 



Multiplying the Hermitian conjugate of Eq. (|E2|) from the right by tpv) gives 
Subtracting Eq. ([E4]) from Eq. (|E3|) gives 



(E4) 



(E5) 



where the Hermiticity of the Dirac matrices has been used. Using the definition of the 
transition current in Eq. ( |63f) , Eq. flE5|) is reveled as the continuity equation 





^ j0» + v ■ jt**) = 



proving the transition-current density to be conserved. 



APPENDIX F: APPLICATION OF GAUSS' THEOREM 



(E6) 



As appendix [E] shows, the transition 4-current, defined in Eq. (0), is a conserved quan- 



tity, 



0. 



(Fl) 



Integrating Eq. ([FT]) over any empty space-time hyper-volume, V, and applying Guass' 
theorem to convert the volume integral into a surface integral over the hyper-surface S 
enclosing V, in general gives, 



daJ^n^ = 



(F2) 
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where the unit 4-vector n M is defined as the outward pointing normal to S. 

For our purposes, it is useful to apply Eq. ( |F2| ) to the space-time region IV, defined in 
Fig. |5| by r± > 0. The closed hyper-surface S enclosing region IV is made of the following 
open hyper-surfaces: (i) t — tf — > +00, (ii) r + = + , r_ > 0, (iii) r_ = + , r + > 0, (iv) 
x — > ±00, and (v) y — > ±00; (see Fig. g). Writing Eq. ( [F2] ) for this surface gives 

p r—00 
0= lim / dr± / J (r,tA 

tf—>OOj J + 00 

-2 f dr ± f° dr„J + {r ± ,T + = + ,r_) 
J J +00 

/r+00 
dr ± j o+ dr + J„(f ± ,r + ,r_ = + ), (F3) 

where use was made of the fact that in any physical situation, i.e. for a square-integrable 
wavepacket, the currents vanish as r — > 00. The hyper- surfaces (iv) and (v) do not contribute 
to the integral of Eq. (|F2j). The factors of 2 arise from the Jacobian relating the original 
differentials to the differentials for the light-front variables, and the negative sign in the 
second and third terms arise because the unit normal vectors h± are directed outside the 
volume V, i.e. J ■ h± = —J±. This completes our proof of Eq. (|55|). 
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FIGURES 

FIG. 1. Schematic diagram depicting a relativistic heavy-ion collision of two charges, and 
Zb-, in the center-of- velocity frame with impact parameter 2b and velocity (3. Lorentz contraction 
is extreme, so the ions are depicted as oblate spheroids. 

FIG. 2. Shown is the scalar component of the Lorentz-gauge interaction, V°, for two different 
energies, (a) 7 = 10 (CERN-SPS energies), and (b) 7 = 100 (RHIC energies), plotted as a function 
of a narrow range of the z-coordinate for t = 0, b = (1, 0), and r±_ = (2, 0). Notice that away from 
the vicinity of z = this interaction is insensitive to the change in the energy of the ion. 

FIG. 3. Same as Fig. ^, except that here the scalar component of the gauge-transformed 
interaction, W°, is plotted. Notice that the gauge-transformed interaction is short-ranged, and 
that the range of the interaction decreases as the ion's energy increases. 

FIG. 4. Schematic diagram of the time- history of a heavy-ion collision in the ultra-relativistic 
limit. Motion in the t — z plane is shown with rj_, and thus b, assumed orthogonal to this plane. 
Ions A and B move toward each other in the z direction with velocity (3 = 1. The dotted lines 
are the projections of the ions trajectories on the t — z plane, which for = 1 coincide with the 
intersection of this plane with the light fronts at z = ±t. 

FIG. 5. The light fronts, i.e. the hypersurfaces defined by t± = 0, divide space-time into 4 
distinct regions: (I) r+ < 0, r_ < 0; (II) r+ > 0, r_ < 0; (III) r + < 0, r_ > 0; (IV) r + > 0, r_ > 0. 

FIG. 6. Depicted is the intersection of the surface S enclosing region IV (shown as the broken 
line) with the r+ — r_ plane. Light-front components of the transition current, J±, are shown 
flowing into region IV at the light fronts, and the time-like component of the transition current, 
Jo, is shown flowing out of IV at the constant, large time tt . 
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FIG. 7. On the left are two terms of our result Eq. (|69|) for the amplitudes, indicated by 
their respective space-time maps: I— > II— ► IV, and I— > III— ► IV. On the right are shown the two 
Feynman diagrams of second-order perturbation theory |[^]with their respective time ordering. 
Our result assumes large 7. Reference [18]assumes small aZ. Exact agreement between the two 
results is obtained in the combined limit. 
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Figure 4, Segev and Wells 




Figure 5, Segev and Wells 
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Figure 7, Segev and Wells 




